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Abstract. We present the results of a systematic search 
for old isolated neutron stars (ONSs) in the direction of 
, two giant molecular clouds in Cygnus (Rift and OB7). 
From theoretical calculations, we expect the detection of 
a large number of ONSs with the PSPC on board ROSAT. 
By analyzing the PSPC pointings in the direction of the 
clouds, we find four sources characterized by count rates 
(~ 10 -3 c s _1 ) and spectral properties consistent with the 
hypothesis that the X-ray radiation is produced by ONSs 
and also characterized by the absence of any measurable 
optical counterpart within their error circle in the digitized 
' red plates of the Palomar All Sky Survey. The importance 
, of follow-up deep observations in the direction of these 
ONS candidates is discussed. The observational and theo- 
retical approach presented here could be fruitfully applied 
also to the systematic search for ONSs in other regions of 
the Galaxy. 

Key words: neutron: stars - X-rays: stars - ISM: indi- 
vidual: Cyg OB7 & Cyg Rift 



1. Introduction 

In the first phases following their formation, neutron stars 
are probably strongly magnetized and can emit intense 
electromagnetic dipole radiation at the expense of their 
rotational energy. These energy losses will produce a non- 
negligible radiation pressure on the surrounding medium, 
inhibiting any possible accretion of interstellar material 
(see e.g. Lipunov 1992). After the radiation pressure has 
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dropped, the flow penetrates inside the accretion radius, 
proceeding unaffected until the magnetospheric (Alfven) 
radius is reached, where the magnetic pressure becomes 
equal to the ram pressure of the infalling gas. At this point 
accretion can continue if the centrifugal acceleration ex- 
erted on the matter flowing along the field lines is smaller 
than the gravitational acceleration: in order to have accre- 
tion from the Interstellar Medium (ISM) a relic magnetic 
field B ^ 10 9 G and a rotational period greater than a few 
seconds are needed (Illarionov & Sunyaev 1975), which 
are not implausible for Old Neutron Stars (ONSs). If this 
condition is satisfied, ONSs should appear as weak, soft 
X-ray sources, as firstly suggested by Ostriker, Rees & 
Silk (1970). 

The detection of ONSs is a long-sought goal of X-ray 
astronomy such that it was included as a possible tar- 
get for the Einstein mission (Helfand, Chanan & Novick 
1980), but the first systematic study of the observability 
of these sources with ROSAT was carried out by Treves & 
Colpi (1991), who found that, in the most favourable case 
of polar cap accretion, thousands of ONSs should appear 
in the ROSAT PSPC All Sky Survey. A complete analysis 
by Blaes & Madau (1993) essentially confirmed the results 
of Treves & Colpi (1991). Zane et al. (1995a) reconsidered 
the detectability of ONSs with ROSAT and found that, 
for polar cap accretion, ONSs may contribute up to 10% 
of the unresolved X-ray background, although in this case 
about 10 sources deg -2 should be observable by ROSAT 
at the sensitivity limit of 10 -3 c s _1 . Finally, in two recent 
investigations Shemi (1995) and Zane et al. (1995b) have 
studied the detectability of nearby ONSs. To date, good 
evidence for ONS candidates has been presented by Stocke 
et al. (1995) and Walter, Wolk & Neuhauser (1996), who 
found that two unidentified X-ray sources detected in ei- 
ther the Einstein Extended Medium Sensitivity Survey or 
in the ROSAT All Sky Survey may actually be powered 
by an accreting ONS. 



2 



T. Belloni et al.: Search for Old Neutron Stars in molecular clouds: Cygnus Rift and Cygnus OB7 



A different approach for the search of ONSs was sug- 
gested by Blaes & Madau (1993) and Colpi, Campana 
& Treves (1993), who indicate as a very promising envi- 
ronment for their detection nearby giant molecular clouds, 
where the density of the ISM is sufficiently high for a single 
neutron star to emit a considerable amount of radiation. 
However, in this case also the absorption of the ISM will 
be enhanced and the detectability will crucially rely on the 
delicate balance between the emitted luminosity and the 
absorbed flux. As shown by Blaes & Madau (1993), Colpi, 
Campana & Treves (1993) and Zane et al. (1995a), some 
molecular clouds in the vicinity of the Sun are expected 
to contain a large number of ONSs emitting at relatively 
high rates and then representing the most favourable sites 
for the observability of these sources. 

On the wake of these studies, we have performed a 
systematic investigation of the X-ray sources detected by 
ROSAT PSPC in the direction of two molecular clouds 
(Cygnus Rift and Cygnus OB7), which should be partic- 
ularly favourable as far as the expected number of ONSs 
is concerned. As we will see, the definition of a systematic 
searching strategy is very important when dealing with 
large number of sources and, to this end, the theoretical 
investigation of the expected emission properties turns out 
to be certainly a key ingredient. 

Section 2 is devoted to the theoretical analysis of the 
number of ONSs detectable by ROSAT PSPC in the di- 
rection of these two molecular clouds and their emission 
properties, whereas in section 3 we will present the search- 
ing procedure and the four ONS candidates which have 
been found in this way. Finally, section 4 contains a dis- 
cussion of these results in connection with the detectabil- 
ity issue of ONSs and a suggestion for possible new deep 
observations to shed further light on the nature of these 
objects. 

2. Theoretical expectations 

2.1. Assumptions 

For a star moving supersonically with velocity v relative 
to the interstellar gas, the amount of mass which is ac- 
creted per unit time is M ~ irr 2 gc nm p v, where r gc = 
2xGM*/v 2 is the gravitational capture radius, M* is the 
mass of the neutron star, n the number density of the ISM 
(n = 1.36 riff, where nn is the hydrogen number density, 
for standard chemical composition) and \ is parameter of 
order unity which encompasses the theoretical uncertain- 
ties on the value of r gc . In the following we take x = 1 ( see 
e.g. Novikov & Thorne 1973). At the typical values of the 
ISM density (outside molecular clouds), it is not clear that 
the fluid approximation applies at the gravitational cap- 
ture radius (see Koester 1976) and the accretion rate could 
be greatly reduced. However, low magnetic fields, frozen 
in the interstellar medium, may effectively couple the gas 
particles. In this case, if the Larmor radius is smaller than 



the gravitational capture radius, the fluid approximation 
can be still considered valid and the accretion rate is rea- 
sonably expressed through the previous formula. The total 
luminosity L emitted by an ONS is then 




where r\ is the relativistic efficiency (77 ~ 0.18 for a neutron 
star radius of ~ 12 km and a mass of ~ 1.4M ) 
and c is the light velocity. If the star moves subsonically 
relative to the ISM, equation (1) remains valid provided 
that the ISM sound speed v s is used in place of v (v s ~ 
10v/(T/10 4 K) kins" 1 ). 

As shown by equation (1), in order to investigate the 
detectability of ONSs, it is crucial to have knowledge of 
the density distribution of the ISM and the velocity distri- 
bution of ONSs in the Galaxy. Calculations of the tempo- 
ral evolution of the ONS distribution function have been 
carried out by many authors (Paczyhski 1990; Hartmann, 
Epstein & Woosley 1990; Blaes & Rajagopal 1991; Blaes 
& Madau 1993; Zane et al. 1995a). Since we are inter- 
ested to investigate the detectability of ONSs in the solar 
neighbourhood and close to the Galactic plane, in the fol- 
lowing we will use the analytic approximation to the ONS 
distribution computed by Zane et al. (1995a). 

The structure of the Local Interstellar Medium has 
been widely investigated: within <~ 50-100 pc from the 
Sun, the gas is tenuous (n ~ 0.05 — 0.07 cm" 3 ) and warm 
(T ~ 10 4 K), although the presence of a number of rela- 
tively high density regions (n ~ 1 cm" 3 ) can rise the ob- 
served column density up to ~ 2 x 10 20 cm" 2 along certain 
line of sights, such as those in the direction of Cygnus Rift 
and Cygnus OB7 (Frisch and York 1983; Paresce 1984; 
Welsh et al. 1994; Diamond et al. 1995). On much larger 
scales, according to Dickey & Lockman (1990), the aver- 
age density in the Galactic plane is n ~ 1 cm -3 with a 
scale height variable with the distance from the Sun, al- 
though there are a number of tenuous (n ~ 0.01 cm -3 ) hot 
(T ~ 10 5 -10 6 K) bubbles. Since Cygnus Rift and Cygnus 
OB7 are very close to the Galactic plane (—5° < b < 5°), 
in the following the ISM will be approximately described 
as a piecewise constant density medium with rin — 0.645 
cm" 3 for d < 100 pc and uh = 1 cm" 3 elsewhere. For 
the sources embedded in the clouds and the background 
sources, the contribution of the cloud material to the total 
absorption will be added. 

A thorough investigation of the emitted spectrum 
turns out to be a further key ingredient in order to cor- 
rectly address the issue of detecting ONSs, since the choice 
of the energy bands where to look for and the absorption of 
the interstellar medium are strongly related to the ONSs 
emission properties. As shown by Alme & Wilson (1973), 
if binary Coulomb collisions between the infalling ions and 
the atmospheric electrons dominate, as it is expected at 
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Table 1. Characteristic parameters for Cygnus Rift and Cygnus OB7 

Cloud volume V c apparent surface S c inner radius n„ outer radius r ou t density n c 

(pc 3 ) (de^) (E£) (E£) (cm- 3 ) 

Cygnus Rift 1.26 xlO b fL2 662 738 29 

Cygnus OB7 l.lOxlO 6 46 739 861 29 



very low accretion rates, the accretion flow can be stopped 
at several Thomson depths and the resulting spectrum can 
be thermalized at a temperature approximately equal to 
the star effective temperature 




where f s is the fraction of the surface area which under- 
goes accretion (see below) . Then, ONSs accreting directly 
from the ISM emit typically in the ultraviolet and soft X- 
ray bands. The low bolometric luminosity (L ~ 10 30 — 10 31 
erg s _1 ) and the softness of the spectrum [T e ff ~ 90 cV) 
explain the difficulty in observing an isolated ONS. Previ- 
ous investigations were restricted to the calculation of the 
emitted spectrum at high luminosities. A detailed numer- 
ical analysis of the spectral properties of unmagnctized 
neutron stars accreting well below the Eddington limit 
has been recently presented by Zampieri et al. (1995). The 
emergent spectrum turns out to be significantly hardened 
with respect to a black body at the star effective tem- 
perature, with a broad maximum shifted toward higher 
frequencies by a factor ~ 3 at L ~ 10 30 erg s~ x . This 
fact is due entirely to the frequency dependence of the 
free-free opacity, for which higher energy photons decou- 
ple at larger depths and temperatures in the neutron star 
atmosphere. 

In the following, we will consider two possibilities for 
the emitted flux: first, black body emission at the neutron 
star effective temperature (see equation 2) which has been 
frequently used in the previous investigations and can be 
regarded as a useful approximation; second, the spectra 
computed by Zampieri et al. (1995). In addition, we as- 
sume that a relic magnetic field is present and that the 
accretion flow is channeled into the polar caps. If we ne- 
glect all the radiative effects produced by the magnetic 
field (on this regard sec e.g. Miller 1992; Shibanov et al. 
1992; Nelson et al. 1995), the main consequence is to limit 
the size of the emitting region by a factor f s = A c /2irr1, 
where A c — 7rr 3 Jta is the area of the polar cap and ta 
is the Alfven radius. For B ~ 10 9 G, as we consider here, 
f s ~ 10 ~ 3 . This fact will produce a hardening of the 
spectrum with respect to the unmagnetized case with the 



same luminosity, since the flux emitted per unit surface is 
Fmag = L/2A C > F unmag = L / Airr 2 , . 

2.2. Expected number of ONSs 

The technique which has been used to calculate the ex- 
pected number of ONSs in the molecular clouds is the 
generalization of a similar procedure introduced by Zanc 
et al. (1995b) and is described in detail in the Appendix. 
With this technique it is possible to calculate also the 
expected number of foreground and background ONSs ac- 
creting from the average ISM, which are seen in the direc- 
tion of the clouds but are not embedded within them. This 
aspect deserves particular notice since, as shall see, they 
are expected to be quite numerous at low sensitivity lim- 
its and could be distinguished from the ONSs embedded 
in the molecular clouds only by means of their observed 
spectral properties. 

The characteristic parameters of the two clouds are 
quoted in Table 1. Results of the expected number of 
sources observable in Cygnus Rift, Cygnus OB 7 and also 
in other spatial regions in the cloud directions are pre- 
sented in Tables 2 and 3 (fifth column) for different values 
of the sensitivity limit S (N^ s = 10 9 , B = 10 9 G). 
These numbers are calculated assuming the spectrum of 
Zampieri et al. (1995) and a black body spectrum (values 
in brackets). The fraction of clouds which has been cov- 
ered by the ROSAT pointings at a given sensitivity de- 
pends mainly on the exposure time. In the second column 
of Tables 2 and 3 we quote the fractional coverage of the 
cloud areas / as a function of threshold. As can be seen, 
the covering increases from ~ 1% at 10~ 3 c s -1 to ~ 5% 
at 1CP 2 c s _1 . To compare the theoretical estimates with 
the actual number of unidentified X-ray sources detected 
in the ROSAT pointings of our sample, for each interval 
of distance we have computed the expected number of de- 
tectable sources according to the following formula 

N(> Si) [dj . d]+l] = (No^sy [dj>dj+l] fi 

+ (Nh h ^)\t d3+l] {h + i-h) + 
+... + 

+ (- /V OJVS r )[^d ;) . + 1 ](/imax - /*max-l)(3) 

where {N h ^g) % ^ d d j is the expected number of ONSs 
above threshold Si detectable in the direction of the whole 
cloud areas and fi is the fractional coverage of the point- 
ings at the same threshold (see Tables 2 and 3). Here 
Si = 1CT 3 c s-\ S 2 = 2 x 1CT 3 c s" 1 , S 3 = 4 x KT 3 
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Table 2. Expected number of ONSs detectable in the direction of Cygnus Rift as a function of threshold (Nqns = W 9 , B = 10 9 
G). Values in the second column (hardness ratio HR) and in the last two columns refer to the Zampieri's spectrum and to a 
black body (values in brackets), dj -h dj+i is the interval of distance considered, (Noffg)^ d ] the expected number of ONSs 
above threshold S detectable in the direction of the whole cloud area (see the Appendix) and N(> S)[ djtdj+1 j the expected 
number of sources in the fraction of the cloud / covered by ROSAT pointings (see equation (3)). 

thresholds' coverage/ hardness ratio H R dj -h dj+i (No l Ns r )[d j ,d j+1 ] N(> S)[ dj ,d j+1 ] 

(cs- 1 ) (%) (pc) 

0.001 0.63 



0.002 2.52 



0.004 3.15 



0.008 4.41 



0.01 5.04 



-0.86-h0.85 (-0.98-^0.69) 


30^100 


3(3) 


0(0) 


0.06-h0.94 (-0.85^0.90) 


100-h300 


64 (43) 


1 (1) 


0.82-h0.97 (0.45-h0.95) 


300-h662 


270 (113) 


7(2) 


0.94-h0.99 (0.87-h0.99) 


662-h738 [Rift] 


297 (215) 


11 (6) 


u.yy ^u.yo— u.yyj 


vqq • op;nn i'vqq • orient 
( oo— zoUU {ioo—Z\Jou) 


OOU \ YLl) 


5 (2) 


-0.81-h0.85 (-0.98-h0.69) 


30^100 


3(3) 


0(0) 


0.07-h0.94 (-0.83^0.90) 


100-h300 


52 (33) 


1 (1) 


0.84-h0.97 (0.53-h0.95) 


300-h662 


167 (73) 


6(2) 


0.95-h0.99 (0.89-h0.99) 


662-h738 [Rift] 


250 (170) 


10 (6) 


n oo (c\ oq ■ n oo\ 
u.yy (^u.yo— u.yyj 


/ oo— zUUU {too— IDoUJ 


loz yoZ) 


4 (1) 


-0.74^0.85 (-0.97-h0.69) 


30^100 


3(3) 


0(0) 


0.17-h0.94 (-0.79^0.90) 


100-h300 


39 (24) 


1 (1) 


0.86-h0.97 (0.62-h0.95) 


300-h662 


96 (46) 


4(1) 


0.95-h0.99 (0.90-h0.99) 


662-h738 [Rift] 


191 (124) 


9(5) 


0.99 (0.98-h0.99) 


738-hl500 (738-hl260) 


47 (17) 


1(1) 


-0.64-h0.85 (-0.96-h0.69) 


30^100 


3(3) 


0(0) 


0.334-0.94 (-0.68^0.90) 


100-h300 


26 (16) 


1 (1) 


0.88-h0.97 (0.69-h0.95) 


300-h662 


51 (27) 


2(1) 


0.96-h0.99 (0.92-h0.99) 


662-h738 [Rift] 


129 (85) 


7(4) 


0.99 (0.98-h0.99) 


738-hllOO (738-h960) 


12 (4) 


1 (0) 


-0.60-h0.85 (-0.96-h0.69) 


30^100 


2 (2) 


0(0) 


0.35-h0.94 (-0.64^0.90) 


100-h300 


23 (14) 


1 (1) 


0.89-h0.97 (0.71-^0.95) 


300-h662 


42 (23) 


2(1) 


0.96-h0.99 (0.92-h0.99) 


662-h738 [Rift] 


110 (74) 


6(4) 


0.99 (0.98-h0.99) 


738-hl030 (738-h880) 


8(2) 


0(0) 



c s-\ S 4 = 8 x 1(T 3 c s-\ S imax = S 5 = 10~ 2 c s" 1 . 
The first term on the right hand side of equation (3) gives 
the number of sources observable in the pointings with 
limiting sensitivity Si, while the other terms account for 
the contributions from pointings at lower sensitivity. Fi- 
nally, the total expected number of sources N(> Si) is 
obtained integrating N(> Sj)^ <j ] over distance and is 
shown in Figure 2 (summing the contributions from the 
two clouds). 

In the fraction of the clouds actually observed by 
ROSAT, we expect the detection of 10-24 sources (the 
exact value depends on the assumption on the emitted 
spectrum) at the sensitivity limit of 1 — 2 x 10~ 3 c s _1 
in the direction of Cygnus Rift and of 5-13 sources in 
the direction of Cygnus OB7. Among these sources, 6-11 
sources are expected to be really embedded in Cygnus Rift 
and 3-6 in Cygnus OB7. The other objects are foreground 
or background sources; in particular, we estimate that 3-8 



and 2-6 foreground ONSs should be detectable in the di- 
rection of Cygnus Rift and Cygnus OB7, respectively. We 
note that at high thresholds 10~ 2 c s _1 the detection of a 
significant number of sources would be expected. 

Then, in agreement with previous investigations (Blaes 
& Madau 1993; Colpi, Campana & Treves 1993; Zane et 
al. 1995a) Cygnus Rift and Cygnus OB7 are expected to 
be particularly favourable sites for the observability of 
ONSs with ROSAT PSPC. In addition, we have found 
that in their direction a significant number of foreground 
and background ONSs should be observable. 

2.3. Expected emission properties 

In principle we could compare directly the theoretical 
spectral distribution with the observed one but, as we 
shall see, the sources selected in our sample have too few 
photons to extract a spectrum. Then, in order to com- 
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Table 3. Expected number of ONSs detectable in the direction of Cygnus OB7 as a function of threshold (Nqns = lo9 > B = 1q9 
G). Values in the second column (hardness ratio HR) and in the last two columns refer to the Zampieri's spectrum and to a 
black body (values in brackets), dj -h dj+i is the interval of distance considered, (No l N°s)ld i ,d i _ l _ 1 ] the expected number of ONSs 



![dj,dj +1 ] 

above threshold S detectable in the direction of the whole cloud area (see the Appendix) and N(> S)[ dj}dj 
number of sources in the fraction of the cloud / covered by ROSAT pointings (see equation (3)). 



the expected 



threshold S 
0.001 



coverage / 
1.55 



0.002 



2.32 



0.004 



3.87 



0.008 



3.87 



0.01 



4.64 



hardness ratio HR 

-0.86-h0.85 (-0.98-h0.69) 
0.06-h0.94 (-0.85^0.90) 
0.82-^0.99 (0.45-^0.99) 
0.95-h0.99 (0.89-^0.99) 
0.99 (0.99) 

-0.81-h0.85 (-0.98-h0.69) 
0.07-h0.94 (-0.83-h0.90) 
0.84-h0.99 (0.53-h0.99) 
0.96-h0.99 (0.90-h0.99) 
0.99 (0.99) 

-0.74-h0.85 (-0.97-h0.69) 
0.17-h0.94 (-0.79^0.90) 
0.86-h0.99 (0.62^-0.99) 
0.96^-0.99 (0.92-h0.99) 
0.99 (0.99) 

-0.64-h0.85 (-0.96-h0.69) 
0.33-h0.94 (-0.68^0.90) 
0.88-h0.99 (0.69-h0.99) 
0.96-h0.99 (0.93-h0.99) 
0.99 

-0.60-h0.85 (-0.96-h0.69) 
0.35^-0.94 (-0.64^0.90) 
0.89-h0.99 (0.71-7-0.99) 
0.97-h0.99 (0.93^-0.99) 



30-hl00 
100-h300 
300-h739 
739-h861 [OB7] 
861-h2100 (861-M700) 

30-hl00 
100-h300 
300-h739 
739-h861 [OB7] 
861^1600 (861-h 1300) 

30-hl00 
100-h300 
300-h739 
739-h861 [OB7] 
861-hl200 (861-h 990) 

30-hl00 
100-h300 
300-h739 
739^861 (739^- 839) [OB7] 
861-h914 

30-hl00 
100-h300 
300-h739 
739-h861 (739-h824) [OB7] 



theor\ 
ONS )[dj 



l] 



3(3) 

64 (43) 
136 (56) 
216 (140) 

67 (22) 

3(3) 
52 (33) 
82 (36) 
167 (105) 

22 (7) 

3(3) 
39 (24) 
46 (23) 
117 (73) 

5(1) 

3(3) 
26 (16) 
25 (13) 
72 (45) 


2(2) 

23 (14) 
20 (11) 

65 (33) 



N(> S) [dj , dj 



pare our theoretical expectations with the observations, 
in each range of distances considered we have simulated 
PSPC spectra by folding the theoretical model with in- 
terstellar absorption, detector effective area and response 
matrix, and calculated the PSPC hardness ratio, defined 

by 

JV(41 - 240)-iV(ll-40) 

A^(41 - 240) + Ar(ll-40) ' 1 1 

where iV(ll - 40) and iV(40 - 240) are the count rates in 
the PSPC channel ranges 11-40 and 41-240, correspond- 
ing roughly to the energy bands 0.1-0.4 keV and 0.4-2.4 
keV. We have repeated the calculation assuming either the 
spectrum computed by Zampieri et al. (1995) or a black 
body at the neutron star effective temperature. Results 
are presented in Tables 2 and 3, where the lowest and 
highest values of the hardness ratio is reported for each 
interval of distance (the highest value refers to the model 
with maximum luminosity). It is interesting to note that, 
for distances above ~ 600 pc for the black body and <~ 300 



pc for the synthetic spectra, the absorption of the inter- 
stellar medium and the cloud material causes the hardness 
ratio to be very close to unity. Then, although ONSs are 
relatively soft sources, above these distances they should 
appear significantly hardened. 

We have estimated also the apparent visual magni- 
tude of ONSs below 100 pc from the Sun and find that 
my <; 32. Then, these sources should lack of optical coun- 
terparts in the digitized red plates of the Palomar All Sky 
Survey (POSS, limiting magnitude m r <~ 20). According 
to recent findings by Blaes, Warren & Madau (1995), for 
ONSs embedded within the clouds and considering polar 
cap emission the reprocessing of the UV-soft X-ray radi- 
ation by the surrounding gas might increase the emitted 
flux at optical wavelengths by 1-2 magnitudes. However, 
even in this case the very low visual magnitude prevents 
any possibility of optical identification for sources beyond 
100 pc. Then, we will use the lack of optical counterparts 
in the digitized POSS plates as a distinguishing criterion 
for selecting ONSs candidates among ROSAT sources. 
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Inserting in equation (A2) the appropriate effective ar- 
eas and frequencies (taken from the world wide web site of 
the Center for Extreme Ultraviolet Astrophysics), we have 
calculated also the expected count rates in the Lex band 
(67-178 A) of the Extreme Ultraviolet Explorer Deep Sur- 
vey (EUVE DS). Below ~ 250-300 pc, the count rate turns 
out to be smaller than 10~ 5 c s _1 even for the more lumi- 
nous sources and no detectable UV flux can be observed 
at Earth. 

Finally, comparing the observed hardness ratios with 
the values quoted in Tables 2 and 3 could be used to dis- 
criminate ONSs among other optically unidentified X-ray 
ROSAT sources. 

3. Data analysis 

3.1. The sample 

We selected from the ROSAT Public Data archive PSPC 
pointings in the direction of the Cyg OB7 and Cyg Rift 
molecular clouds. A description of the satellite and the de- 
tector can be found in Triimper (1983) and Pfeffermann 
et al. (1986). The 28 analyzed areas are shown in Fig- 
ure 1, overlaid onto the CO contour maps from Dame & 
Thaddeus (1985). The relevant data on the pointings are 
reported in Table 4. A few pointings lie just outside the 
contour of the clouds in Figure 1, but have been never- 
theless included in the sample. The total exposure time 
is 126603 seconds. Only the central region of the PSPC 
(radius 20') has been considered for the analysis, to avoid 
problems due to the poor knowledge of the point spread 
functions outside that region. The areas of the two clouds 
covered by ROSAT pointings are 2.1 square degrees (Cyg 
OB7) and 5.6 square degrees (Cyg Rift), which correspond 
to 4.6% and 5.0% of the total areas of the clouds, respec- 
tively (adopting the values quoted in Dame & Thaddeus 
1985). Because of the nature of the pointings, the frac- 
tional coverage of the cloud areas is limited and varies 
with threshold (see Tables 2 and 3). Although at the low- 
est sensitivity limits (1-2 10~ 3 c s _1 ) only 1-2% has been 
covered, the theoretical expectations (see again Tables 2 
and 3) indicate that this could be sufficient for detecting 
ONS candidates. Moreover, from the systematic analysis 
of these data a search methodology can be outlined. 

3.2. The analysis procedure 

Each pointing has been extracted from the ROSAT Pub- 
lic Archive at MPE Garching and analyzed separately by 
running a set of semi-automatic programs. The standard 
detection technique in the EXSAS package (Zimmcrmann 
et al. 1994) has been applied. First a background map is 
produced by removing all possible sources (identified by 
means of a sliding window technique) and running a two 
dimensional spline fit to the data. Then a Maximum Like- 
lihood (ML) algorithm is applied (see Cruddace, Hasinger 
& Schmitt 1988) to detect significant deviations from the 



Table 4. Observation log. The first two columns are the iden- 
tifier of the pointing in this paper and the ROSAT observation 
number. Galactic coordinates are in decimal degrees. Nx is the 
number of detected X-ray sources and Ngsc is the number of 
GSC stars used for the boresight correction (noB indicates one 
star or less, see text) . 



Label 


DAD _LL 

rUJrt 7f 


rjXp. [s) 


i 

in 


Oil 




Ng5C 






Cyg OB7 










A 


200427 


382.0 


91.4 


1.1 


1 


noB 


K 


500143 


8260.0 


93.3 


6.9 


9 


4 


L 


500254 


4297.0 


93.7 


-0.3 


1 


noB 


S 


300033 


4331.0 


89.8 


-0.1 


5 


noB 


T 


300033.1 


2854.0 


89.8 


-0.1 


6 


3 


W 


500195 


9137.0 


94.0 


1.0 


7 


5 


b 


900194 


4166.0 


89.4 


-0.7 


3 


2 






Cyg Rift 










B 


200057 


1846.0 


76.6 


1.4 


4 


3 


C 


200058 


1931.0 


76.6 


1.4 


2 


2 


D 


200062 


14030.0 


76.6 


1.4 


17 


14 


E 


200063 


2073.0 


76.6 


1.4 


3 


3 


F 


200109 


3587.0 


80.2 


0.8 


11 


9 


G 


200114 


3508.0 


85.7 


1.5 


1 


noB 


H 


200564 


8950.0 


74.3 


1.1 


5 


3 


I 


500084 


2470.0 


69.7 


1.0 


2 


2 


J 


500085 


2879.0 


83.0 


-0.3 


2 


2 


M 


900157 


4054.0 


78.2 


2.0 


2 


1 


N 


900158 


4123.0 


76.0 


0.4 


1 


noB 


O 


200154 


1261.0 


75.8 


1.3 







P 


200412 


1265.0 


78.2 


0.1 


1 


noB 


Q 


200412.1 


2228.0 


78.2 


0.1 


2 


noB 


R 


200702 


4409.0 


85.7 


-0.3 


7 


5 


U 


400157 


2996.0 


77.0 


3.2 


3 


2 


V 


500124 


5979.0 


70.7 


1.2 


1 


noB 


X 


500203 


3641.0 


84.2 


-0.8 


1 


noB 


Y 


500206 


10323.0 


78.4 


2.8 


13 


6 


z 


500207 


9319.0 


77.9 


3.5 


8 


5 


a 


500208 


2304.0 


77.5 


4.0 


2 


noB 


Total 




126603 






120 



estimated background distribution. The threshold for de- 
tection has been set to the conventional value of ML=10, 
which corresponds to a chance detection probability of 
4.5 xl0~ 5 for a single trial. This technique is repeated for 
the data in three different PHA channel ranges: total (T: 
channels 11-240, 0.1-2.4 keV), soft (S: channels 11-40, 
0.1-0.4 keV) and hard (H: channels 41-240, 0.4-2.4 keV). 
In case of a detection in more than one band, only the 
band with the highest value of ML has been considered. 

It is known that the ML technique becomes problem- 
atic in presence of large regions of extended emission, as it 
is likely to find in the pointing directions considered here, 
lying on the galactic plane. Therefore, all the detected 
sources have been visually inspected in the X-ray images 
and obvious spurious detections have been discarded. 
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Fig. 1. Map of the analyzed ROSAT PSPC pointings. Each circle is 20' radius. The black circles correspond to pointings 
containing ONS candidates (see text). The contours of the clouds are from Dame & Thaddeus (1985). 



Fig. 2. The log iV-log S distribution for all the sources de- 
tected by ROSAT (continuous line) and all the optically 
unidentified sources detected by ROSAT (dashed line). Also 
shown is the expected number of ONSs, N(> S), in the frac- 
tion of the clouds covered by ROSAT pointings (dash-dotted 
line; Zampieri's spectrum). 

In the cases where two or more pointings were 
coaligned, they were analyzed separately: the sources have 
been cross-correlated and in case of coincidences only the 
source with the highest detection ML has been considered. 
Eleven sources have been eliminated in this way, bringing 
the total number of sources to 109. The logiV-logS* dis- 



tribution for the sources detected by ROSAT is shown in 
Figure 2. There are a few strong sources, but most of them 
are below 0.01 c s _1 . Note that the targets of the pointings 
have not been removed. 

3.3. Optical identification 

The positions of the detected X-ray sources had to be 
corrected for an offset due to a residual systematic uncer- 
tainty in the position determination (the so-called bore- 
sight correction). In order to do this, preliminary identi- 
fications have been performed using the HST Guide Star 
Catalog (GSC; Lasker et al. 1990). The X-ray sources and 
GSC stars have been cross-correlated and the best shifts to 
the X-ray positions have been determined by means of a 
Maximum Likelihood technique. The shift so obtained has 
been considered only if at least two sources were identified 
in this way (with the exception of pointing S, where one of 
the two sources identified with GSC entries had multiple 
bright optical counterparts). For the remaining pointings, 
a 6" (la) error has been quadratically added to the error 
radii of the detected sources (Kiirster & Hasinger 1992). In 
Table 4 a summary of the X-ray detections and boresight 
corrections is reported. An X-ray source is considered to 
be identified if a GSC star falls into its 90% positional er- 
ror circle. In case of more than one optical identification, 
the brightest counterpart has been considered. 

67 X-ray sources were identified with entries in the 
GSC, which has a loosely defined magnitude limit of 
my ~ 14. For the remaining 42 sources, we used the Dig- 
itized Sky Survey (DSS; Postman et al. 1995), which con- 
sists of the POSS red plates. The magnitudes of the bright- 
est stars in the error boxes have been estimated by using 
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the empirical magnitude-diameter relation of King & Raff 
(1977) appropriate for red POSS plates. The values so ob- 
tained have a large uncertainty (up to one magnitude), 
but it is sufficient for our purpose. Only in three cases 
the possible optical counterpart was not bright enough to 
allow an estimate of the magnitude, the image being not 
saturated. In Figure 3 we show the distribution of mag- 
nitudes of the optical counterparts for the two samples, 
GSC and DSS. As one can observe there is a considerable 
number of bright stars, since there are bright OB associ- 
ations in the direction of most of the pointings; indeed a 
large number of the counterparts for which we were able 
to obtain a spectral type were of types O and B. The pres- 
ence of strong optical candidates for many X-ray sources 
ensures the robustness of the boresight process. 

After the optical screening, 7 sources were left for 
which no optical counterpart could be seen in the POSS 
red plates, to a limiting magnitude m r ~ 20. Two of these 
were coincident with the original targets of the pointings 
and were therefore excluded. The brightest of the remain- 
ing sources was subsequently identified with Cyg X-3, so 
4 final candidates remained. 



3.4- ONS candidates 

For the four sources with no optical counterpart, PSPC 
count rates in the different bands have been extracted 
manually from the data, with an accurate estimate of the 
background, in some case contaminated by extended emis- 
sion. The relevant information about the candidates (two 
in each cloud) can be found in Table 5. As it can be seen 
from Table 4, all four sources come from pointings with 
long exposure times and with a large number of GSC iden- 
tifications, which ensures an accurate boresight correction. 
Source OB7-2 lies just outside the outer rim of the cloud 
(see Figure 1), but we decided not to discard it from the 
sample. 

Source Rift-2 might have counterparts at different 
wavelengths. The POSS plate shows a very faint ring- 
like structure (~ 20" diameter) in the X-ray error box, 
reminiscent of a supernova remnant or a planetary neb- 
ula. The position of the source is compatible with that of 
the IRAS source 20179+4102, which has a 100 fim flux 
of 75 Jy (IRAS Catalogs and Atlases 1985). Putting the 
distance at 0.7 kpc, this translates roughly to a luminos- 
ity of ~ 10 35 erg s _1 . At the same distance the radius of 
the optical ring would be <~ 0.1 pc. If the ring and the 
IRAS sources were the same object, it would be too small 
and faint to be a supernova remnant, while the planetary 
nebula hypothesis could not be ruled out. In absence of a 
clear identification we do not exclude the source from our 
candidates list. 



Fig. 3. Distribution of optical magnitudes for the possible 
counterparts of the X-ray sources. GSC magnitudes are in the 
V band, POSS magnitudes in the red (see text). The lack of 
POSS stars at magnitude 19 is a probable artifact of the pro- 
cedure for magnitude estimate. 

4. Discussion 

As can be seen from Tables 2, 3 and 5, our predicted num- 
ber of sources is larger than the number of unidentified ob- 
jects detected by ROSAT. In particular, Figure 2 shows 
the expected number of ONSs and the observed number of 
unidentified objects in our sample as a function of thresh- 
old. By comparing the two curves, it is clearly visible that 
predictions are too "optimistic", although it should be 
stressed that, by purpose, theoretical calculations have 
usually been performed under the most favourable as- 
sumptions. The fact that the observed number of ONS 
candidates is smaller by a factor 3-10 indicates that the 
range of parameters of our model needs to be more con- 
strained. As already discussed, there are several sources of 
uncertainties which may influence theoretical predictions, 
such as the velocity distribution and density of ONSs, the 
strength of the relic magnetic field (related to the possi- 
ble inhibition of accretion), the neutron star birth rate, 
the total number of ONSs, the indetermination in the ex- 
act value of the gravitational capture radius r gc and the 
validity of the fluid approximation for the average ISM. 
Scattering by molecular clouds and spiral arms (dynami- 
cal heating; Wielen 1977) may boost a fraction of low ve- 
locity stars to higher speeds, reducing the number counts 
up to an order of magnitude (Madau & Blaes 1994). Also, 
evidence has been found (Lyne & Lorimer 1994) for a neu- 
tron star mean birth velocity much higher than the one 
obtained by Narayan & Ostriker (1990). Recently, Blaes, 
Warren & Madau (1995) have shown that preheating of 
the ambient gas to temperatures higher than 10 4 K may 
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Table 5. ONS candidates. The columns are: source name, pointing identifier, celestial coordinates, 90% error radius, Maximum 
Likelihood of existence, PSPC total count rate (channels f f-240) and hardness ratio (defined as HR = (H-S)/(H+S), where S 
and H are the counts in the channel intervals lf-40 and 4f-240 respectively. 



Name 


P. 


a (2000) 


5 (2000) 


Ar(") 


ML 


c ksec" 1 


HR 


Rift-f 


D 


20' 1 19 m 27.96 s 


+38°38'27.6" 


10.5 


24.7 


2.7±0.6 


0.29±0.23 


Rift-2 


Y 


20 h 19 m 47.05 s 


+41°12'0f.3" 


19.2 


10.9 


3.2±0.9 


>0.28 


OB7-1 


K 


20 h 53 m 07.35 s 


+55°f3'29.0" 


13.4 


17.7 


1.5±0.5 


>0.44 


OB7-2 


W 


21 h 25 m 31.12 3 


+51°48'31.5" 


9.2 


35.2 


2.4±0.6 


>0.84 



significantly reduce the accretion rate, hence decreasing 
the number of detectable objects. 

Our ONSs searching strategy was aimed to select X- 
ray objects in molecular clouds which, within their error 
boxes, do not have any optical counterpart in the POSS 
plates. Additional information is clearly needed to assess 
the true nature of these objects. The detection of an X- 
ray pulsation in the range of 1-1000 s will clearly indicate 
that we are observing an ONS; however high magnetic 
fields are needed which could prevent accretion and could 
harden substantially the emitted spectrum (see e.g. Nel- 
son et al. 1995), moving the peak emissivity outside the 
ROSAT energy range. In absence of X-ray pulsations, a 
spectral criterion could be used to individuate ONSs: in 
the case of low magnetized or unmagnetized ONSs we ex- 
pect an X-ray spectrum which show a significant harden- 
ing with respect to a black body at the neutron star effec- 
tive temperature, as discussed in Zampieri et al. (1995). 
An additional feature could be a high column density, tes- 
tifying that the X-ray object lies within (or beyond) the 
molecular cloud. As can be seen from Tables 2 and 3, 
the expected hardness ratios for ONSs embedded in the 
clouds or beyond them are near to unity. In this respect, 
the observed HR of three of our candidates are consistent 
with the hypothesis that the sources are located within 
the clouds, whereas one source (Rift-1) cannot certainly 
lie within the cloud because of its measured value of the 
hardness ratio HR = 0.29 ±0.23 (see Table 5). However, it 
should be noted that its spectral properties are consistent 
with those of a foreground ONS located at <: 100-300 pc 
(see also Table 2a, £ = 2 x 10~ 3 c s _1 ). We checked that, 
at this distances, Rift-1 should not be detectable even at 
the lowest sensitivity limits 10~ 3 c s" 1 ) of EUVE DS. 

From the lack of optical counterparts in the digitized 
POSS plates, we can estimate a lower limit for the ob- 
served X-ray to optical flux ratio (Fx/Fy) for our ONS 
candidates. The limiting magnitude in the regions of the 
four candidates is m r ~ 20; using the conversion to magni- 
tude uiy by King and Raff (1977) we can adopt this limit 
as a conservative limit on my, since only blue objects 
would break this assumption and the bluest unabsorbed 
stars have B-V=— 0.3. For the conversion between ROSAT 
countrate and flux we adopted the relation by Fleming et 



al. (1995), appropriate for sources of stellar nature. The 
resulting values for log Fx jFy are between -0.4 and +0.1. 

Stocke et al. (1995) use the value of Fx jFy as a chief 
criterion for determining the nature of the X-ray emitter. 
In addition to compact objects, high Fx jFy are associ- 
ated also with BL Lacs, cluster of galaxies, young super- 
nova remnant and heavily extincted pre-main sequence 
stars. While an extragalactic origin seems unlikely for 
sources on the galactic plane in the direction of molecular 
clouds, the other possibilities still cannot be excluded. 

Follow up deep observations will be necessary to fur- 
ther investigate the nature of these ONS candidates. In 
particular, deep X-ray pointings in the direction of these 
four sources would be very useful to measure with more 
precision the X-ray spectral properties. Also, UV obser- 
vations in the LeX band of EUVE at the lowest sensitivity 
limits would be desirable to verify the absence of ultravi- 
olet emission. Finally, deep optical observations are cer- 
tainly needed to search for the presence of possible fainter 
counterparts within their error circle and to put a better 
constraint on the Fx jFy ratio. 

A. Appendix: expected number of ONSs 

In the calculation of the expected number of ONSs we 
have used the analytic approximation to the ONS distri- 
bution computed by Zane et al. (1995a), which is based 
on the evolution of the F population (characterized by the 
lowest value of the mean velocity) of model b of Narayan 
& Ostriker (1990) and on the assumption of a uniform 
spatial distribution of neutron stars: 



dN ONS 
dVdv 



noNS 



= n Q NS 



dG(v) 



dv 



dv 



1 + (i^kHis- 1 ) 3 - 3 



(Al) 



where dV is the spatial volume integration element, 
n NS = 3 x 10~ 4 (N^g/W 9 ) pc~ 3 is the number den- 
sity of ONSs averaged in the local region around the Sun 
(r < 2 kpc, \z\ < 200 pc), Nq^ s is the present total num- 
ber of ONSs in the Galaxy and G(v) is an analytic fit to 
the cumulative velocity distribution derived by Zane et al. 
(1995a). 
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The count rate (measured at Earth) in a certain spec- 
tral interval (y\, v 2 ) of a star which emits a monochro- 
matic flux of radiation F v is 

^(7) e-°» NB A v du (A2) 

where d is the distance of the source, ov is the absorption 
cross-section of the ISM (Morrison & McCammon 1983) 
and Nh the hydrogen column density. A v is the effective 
area and Av = v 2 — v\ the bandpass of the detector (y\ = 
0.1 keV and v 2 = 2.4 kcV for ROSAT PSPC). 

Clearly at a certain distance dj an ONS will be de- 
tectable if its count rate N is above the sensitivity limit S 
of the detector. Then, at dj there exists a limiting value 
of the star luminosity, {L min )d. which depends on the 
emission properties and the absorption of the ISM, be- 
low which an accreting neutron star does not give rise 
to a count rate above the threshold of the detector and 
hence is not observable. From equation (1) this translates 
directly into an upper limit for the star velocity, (y max )d. . 

Alternatively, a star with a given luminosity L, or ve- 
locity v, will be observable up to a maximum distance 
d m ax{v) at which N goes below S. Then, the total num- 
ber of ONSs which can be observed in a certain inter- 
val of distance [dj, rfj+i] and within a solid angle £ 
can be calculated from equation (Al) summing up all 
the neutron stars that are contained within the volume 
V(v) = S[rf^ a;E (w) — dj]/3 and integrating over v. 

(Afe r )[^ +1 ] = n ONS / V(v) — dv + 

Jv j+1 av 

+n NsV(vj +1 )G(v j+ i) (A3) 

In equation (A3) the integral goes from the maximum ve- 
locity of detectability at distance dj+i, Vj +1 = (wmm)^, 
to the maximum velocity at dj, Vj = (tWa;)^ , whereas 
the second term on the right hand side accounts for all 
the stars with luminosity above threshold throughout all 
of the spatial volume considered. The integral has been 
evaluated numerically using a Lobatto quadrature and 
dmax(v) has been interpolated at the appropriate value 
of v from a table of entries previously calculated (an ex- 
tensive use has been done of a local cubic interpolation 
procedure). We note that for each threshold there exists 
an absolute upper limit to the distance of detectability 
which corresponds to the star with the minimum accre- 
tion velocity v — v s and maximum emitted luminosity. 
We note also that in each interval of distance [dj, 
the ISM density must be constant in order to ensure the 
regularity of the function d max (v) and hence the correct 
computation of the integral. 

We have estimated the spatial boundaries rj„ = r c — 
Ar c and r out = r c + Ar c of Cygnus Rift and Cygnus OB7 
(r c is the position of the cloud center) by inverting for the 
cloud width Ar c the expression V c = £ c (r 3 ^ — rf n )/3, 



where V c is the cloud volume and S c its apparent angular 
surface (taken from Dame et al. 1987; see also Table 1). 
Using these boundaries, the expected number of ONSs in 
each cloud has been evaluated from equation (A3). More- 
over, with the same technique, it is possible to calculate 
also the expected number of foreground and background 
ONSs accreting from the average ISM, which are seen in 
the direction of the clouds but are not embedded within 
them. 
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